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Abstract 

Long sequences of successive direct (projective) measurements or observations of a few 
“uninteresting” physical quantities of a quantum system may reveal indirect, but precise and 
unambiguous information on the values of some very “interesting” observables of the system. In 
this paper, the mathematics underlying this claim is developed; i.e., we attempt to contribute 
to a mathematical theory of indirect and, in particular, non-demolition measurements in 
quantum mechanics. Our attempt leads us to make novel uses of classical notions and results of 
probability theory, such as the “algebra of functions measurable at infinity”, the Central Limit 
Theorem, results concerning relative entropy and its role in the theory of large deviations, 
etc. 


1 Introduction: Purpose and scope of paper 

This paper is devoted to a study of the theory of indirect measurements and observations in 
quantum mechanics, m- Our main aim is to develop a general perspective on questions that, 
in very general terms, can be formulated as follows: What sort of information on a quantum 
system, S, can be extracted from long, time-ordered sequences of direct (so-called projective) 
measurements or observations of a single physical quantity of S represented by a self-adjoint 
operator, V; (or of a finite number of physical quantities, A' := {AJ,..., X r }, of S)? What do 
such sequences of direct measurements/observations reveal about facts concerning S ? What does 
quantum mechanics tell us about the time evolution of quantum systems subjected to repeated 
measurements/observations? These are fairly fundamental questions about quantum theory that 
call for general answers arrived at in an exploration of general concepts and structure that goes 
beyond the analysis of special examples. The analysis presented in this paper is inspired, to some 
extent, by the work of experimental groups, see, in particular m, and by various theoretical 
papers, see EH, ML 0, 0, 0, and references given there. Although our main interest concerns 
conceptual aspects of quantum mechanics, the main effort that has gone into this paper concerns 
the study of simple models (see, e.g., [3.) describing a concrete physical situation. This has the ad¬ 
vantage that it renders our paper comprehensible for readers who are not used to abstraction. Yet, 
our analysis is intended to illustrate some - we believe novel - general insights: The mathematics 
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underlying the theory of long, time-ordered sequences of direct measurements of some physical 
quantities of a quantum system S turns out to be closely related to the one underlying the theory 
of equilibrium (Gibbs) states of long chains of classical lattice gases (or classical spins), with time 
of the quantum system corresponding to ID physical space of the lattice gas and sequences of out¬ 
comes of direct quantum-mechanical measurements corresponding to configurations of particles 
in the lattice gas. Results of statistical mechanics, such as the equivalence of different ensembles 
in the thermodynamic limit, or the disjointness of equilibrium states corresponding to different 
values of some thermodynamic parameters, such as the temperature or a chemical potential, have 
close cousins in the quantum theory of repeated measurements. For example, the observation of 
an event of finite duration in a quantum system corresponds to the appearance of short-range 
order in some bounded spatial region of the lattice gas, and an “eternal property” of a quantum 
system (i.e., a property of S observed in a non-demolition measurement) corresponds to a specific 
value of an order parameter labeling different equilibrium states of the lattice gas. Mathematical 
methods, such as large-deviation theory, concentration-of-measure estimates, hypothesis testing, 
etc., can be transferred from statistical mechanics to quantum theory. 

In this paper, we will study a quantum system, S , with the following properties: 

1. For all practical purposes, S can be considered to be an isolated (“closed”) system. This 
means that, in the Heisenberg picture, the equations of motion of operators representing 
physical quantities of S take the form of Heisenberg equations of motion in which the “Hamil¬ 
tonian” of S appears. 

2. S' is the composition of two subsystems, P and E, where P is the system we actually wish 
to study, while E consists of all the experimental equipment - probes, detectors, and other 
measuring devices - used to observe P. Clearly, P and E interact with one another; so neither 
P, nor E are isolated systems; while S = P\IE is isolated if E is chosen appropriately large; 
see e.g. |T2] . 

We propose to analyze a concrete model system where P is a quantum dot in a semi-conductor 
device containing a component, P, close to E that can bind up to N < oo electrons. A sketch of 
this system is given in Figure 1. 



Electrons can enter into the component P or tunnel out of it and into a component P' , (i.e., P 
can be ionized). However, processes changing the number, v, of electrons bound by P are very slow 
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(in a sense to be specified below). The experimental equipment E consists of a conducting channel 
close to the region where P is located, with electron detectors, Dr and Dr, at both ends of the 
channel, and of an electron gun that shoots electrons into the channel at regular time intervals; 
(e.g., one electron every r seconds, where r is such that, typically, there is only one electron in 
the channel). The electrons in the channel connecting Dr and Dr serve as probes to observe P 
indirectly: the rates at which electrons in the channel hit Dr or Dr, respectively, depend on the 
number, v, of electrons bound by P. These latter electrons give rise to a “Coulomb blockade” in 
the branch of the channel leading towards Dr. the larger v, the more likely it is that an electron 
propagating through the channel will hit the detector Dr. 

The only direct observation of S that can be made is to see whether an electron in the channel 
triggers a click of either Dr or Dr. In an idealized (simplified) mathematical description of S, 
this observation will be represented by an operator, X, with 

x =%®(; °i) E - < li > 

with the property that X has the eigenvalue 1 on all channel states describing only one electron 
propagating through the channel that will hit Dr and the eigenvalue —1 on all channel states 
describing one electron hitting Dr. Failed measurements will not be recorded. When an electron 
is scattered onto a detector it disappears for ever, and the original state of the detector is restored 
after a time much shorter than the cycle time r of the electron gun. The physical quantity 
represented by the operator X is measured around times r, 2r, 3r,.... These measurements are 
projective measurements, as described in El- If the electron gun and the detectors function 
perfectly, only the eigenvalues £ = ±1 are observed in every measurement. We thus denote by 
= ±1 the outcome of the j th successful measurement of X; (with = 1 if Dr clicks, and 
= — 1 if Dr clicks). As is evident from Figure 1, the probability of the outcome = 1 in the 
j th measurement of X depends on the number v of electrons bound by the component P of the 
quantum dot P close to the conducting channel. Because of the effect of the Coulomb blockade 
on the motion of electrons in the channel, the probability of the outcome = 1 increases with 
increasing v. It is convenient to introduce an operator 


<yF := c/tp ® Ip/ (S> Ip (1.2) 

whose eigenvalues v £ {0,1, 2,..., N} correspond to the number of electrons bound by the com¬ 
ponent P. Clearly, at a fixed time, ,/F commutes with X. We assume that the Heisenberg time 
evolution of ,zF is very slow, in the sense that the norm of 

(i-3) 

is tiny for all t, where ../F ( t ) is the Heisenberg-picture operator corresponding to ^Y at time t. 
We begin by studying the special case where ./F commutes with the time evolution, i.e., ,/F(f) is 
independent of t. One then speaks of a “non-demolition measurement” of the number of electrons 
in P. 

For our model, a measurement protocol of length k consists of a sequence 

i (fc) = (i-4) 

k = 1,2,3,...., of outcomes, = ±1, of k successful successive projective measurements of 
X; see, e.g., |T5]. (Recall that unsuccessful measurements of A' are not recorded, and we will 
henceforth ignore them.) Quantum mechanics enables us to calculate the probability, denoted by 
A t w(£ii£ 2 >—>£*)> of a measurement protocol ^ k \ given the state oj of the system S just before 
measurements of X have started. 

The question we wish to address is what a measurement protocol tells us about the number, 
v, of electrons bound by P. The answer is actually fairly obvious. Given v, one calculates the 
repulsive potential acting on an electron in the right branch of the channel near the region where P 
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is located. Simple quantum-mechanical calculations then yield the Born probabilities, p(L\v) and 
p(R\is), for an electron in the channel to be scattered onto Dr or onto Dr, respectively. Clearly, 

p{L\v) + p(R\v) = 1, (1.5) 


and if the detectors function perfectly either Dr or Dr will click, each time the electron gun has 
shot an electron into the channel. Given a measurement protocol let 

fi L k \0 ■■= -jTZ-fij e {/ + 1, k} I = 1} (1.6) 

be the frequency that an electron - out of k — l electrons shot into the channel after the I th shot 
of the gun - is scattered onto Dr, and let fR ~ k \£) be the frequency that an electron is scattered 
onto Dr, with | := (£i,£ 2 , ■••) = (&)»eN- Then 


f { L’ k \0 + fR k \0 = 1 - ( 1 - 7 ) 

The law of large numbers suggests that the limit 



exists for any choice of l, a.s.. The value of /{°^(0 is independent of l and of the first m 
measurement outcomes, £i,...,£ m , for any m < oo. Such a function is called an “observable at 
infinity”, or a “pointer observable”. Now, if the state of P happens to be an eigenstate of the 
operator JV corresponding to the eigenvalue v (i.e., v electrons are bound by P, and tunneling 
between P and P' = P\P is suppressed) then one expects that 

f { L %(i)=p(L/R\v), (1-9) 


as quantum theory predicts. Once we understand how fast the limit k —> oo on the right side of 
(1.8) is approached, and assuming that p(L\v) separates points of the set {1, ...,7V}, meaning that 


min \p(L\vi) — p(L\v 2 )\ =: k > 0, (1.10) 


then we will have found the answer to the question what a measurement protocol is likely to 
tell us about the number v of electrons bound by P, assuming that k is large enough: v can be 
inferred from (£), (with l fixed), by using (1.8)-(1.101, with an error rate that tends to 0, as 
k —> oo! 

Actually, things are a little more complicated, but only slightly. If P consists of two or more 
spatially separated components joined by tunneling junctions, with only P close to the channel 
connecting Dr and Dr, then the quantum state of P is usually not an eigenstate of the operator 
jV even if that state is pure. The probability for jV to have the value v £ (1,2,..., N} - when 
measured - is determined by the quantum state of P according to Born’s rule. A rather straight¬ 
forward argument shows that, after the passage of a large number of electrons through the channel 
joining Dr and Dr, the state of P approaches an incoherent mixture of eigenstates of , i.e., of 
states with a definite number of electrons bound by P. This phenomenon is called “ decoherence”. 
Once decoherence has set in, the claim after (1.101 becomes valid. 

Our assumption that Yf commutes with time evolution, i.e., that a non-demolition measure¬ 
ment of is carried out, is only justified in the limit where the probability for an electron to 
tunnel into or out of P tends to 0. While this tunneling probability may be very small, it is usually 
non-zero. We should therefore ask what can be said about the time evolution of the eigenvalue 
v(t) of the operator ^ Y(t ) (with t = time ) inferred from long but finite measurement protocols 
recording the outcome of repeated direct measurements of the observable represented by the op¬ 
erator X. We will show that the average number of electrons bound by P during time intervals 
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of length T»r, as revealed by very long protocols of measurements of X, can be described by 
some stochastic jump process, (v(t)) t= jTj= 0 , 1 , 2 ,...) on the spectrum {1, ...,N} of the operator jV, 
where T is determined by the speed of convergence in (1.8 1 and the amount by which ^V(t) varies 
over the cycle time r of the electron gun, with T —> 00 in the limit where ^V{t) is independent of t 
(non-demolition measurement of ^V). A fairly elementary result on the process 0 , 1 , 2 ,... is 

presented in Section[4j More refined results on the properties of the jump process {v(jT))j= 0 , 1 , 2 ,... 
will appear in forthcoming work. 

We pause to comment on the analogy between the quantum mechanics of sequences of projec¬ 
tive measurements, as described above, and the statistical mechanics of classical lattice gases: A 
measurement protocol corresponds to a configuration of particles of the lattice gas; (£j = lf> 
a particle occupies site j £ Z+, = — 1 site j is empty). The quantum-mechanical prob¬ 
ability, Mw(£i>£2,-->£fc)> °f observing the sequence = (£ 1 ,.••,£&) of measurement outcomes 
corresponds to the probability of the corresponding configuration of particles, as predicted by the 
Gibbs equilibrium measure of the lattice gas. The frequency (£) corresponds to the density 

of particles in the subset {1,2, ...,fc} of the one-dimensional lattice Z + , and corresponds 

to the infinite-volume density of particles. The number v of electrons bound by P corresponds 
to the chemical potential of particles in the lattice gas. The equivalence between the canonical 
and the grand-canonical ensemble in the thermodynamic limit of the lattice gas corresponds to 
the statement that has a sharp value, p(L\v), and that the fluctuations of /|°’^(£) are °f 

order 0(k~ s), as k —► 00 ; etc. 

We conclude this introduction by adding a little more precision to the statements sketched 
above. Let H denote the space of infinitely long measurement protocols (equipped with the a- 
algebra generated by the cylinder sets 


'■= 1C e s I Zji = Vj t ,i = 1 (1.11) 

for arbitrary choices of (r]j 1 ,..., m < 00 ). Given £ £ S, we denote by the sequence 

of measurement outcomes (£*+i, £;+ 2 j £fc), l < k, and, apparently, ^ 0,fc) . As has been 

indicated above, the choice of a state, ui, of the quantum system S = P V E equips the measure 
space 5 with a probability measure, p u , defined on the cr-algebra generated by the cylinder sets 
■■■■ rij m ). It associates with a protocol the generalized Born probability 

( 1 - 12 ) 


for the first k projective measurements of the observable X (carried out at times ~ jr, j = 
1,..., k) to yield the values £ 1 ,..., The unique expression for /u w (£i,..., £*,) provided by quantum 
mechanics was first found by Schwinger in [26] and rediscovered by Wigner m and, most certainly, 
by very many other theorists. It is recalled in Sect. 2.1 Eq. (2.4). For consistency of // w , it is 
necessary that 

'y ^ /£i;(£lj •••) Cfc) £fc+l) = ^fc)> (1.13) 

ik+1 


with Huj(0 = lj so that p u is a probability measure on S. These conditions hold always true 
in quantum mechanics. One then observes that 




VjJ) 


^ ■■■! €jm)- 



(1.14) 


The general theory of direct measurements/observations in quantum mechanics requires that 
Hu, also satisfies a certain “decoherence condition”, which says, roughly speaking, that subsequent 
measurements of the observable X are “independent” of one another, (which does not mean that 
the measurements outcomes £ 1 , are uncorrelated). For the purposes of this introduction 
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we do not need to know precisely what is meant by decoherence and why it is important; (but see 
Sect. [ 3 ] and, c.g., mm)- A consequence of decoherence is that 


^ ) Pui(£,lj £fc+l> ) £,k+r) — pui{^k+ 1) £fc+r)j 

Ci, 


(1.15) 


(possibly only up to a tiny error), for arbitrary k,r £ N, where the right side of (1.15) is the 
generalized Born probability for the outcomes £fc+i, £k+r in measurements k + 1 + r, 

assuming that measurements 1,..., k have not been carried out; (there have been no gun shots at 
times r, 2r,..., hr). 

Next, we define certain subsets of 5: 


B^GO := {£ e 3 I \f { L k) (0-p{L\v)\ < s k -i}, l < k, 


(1.16) 


where, we recall, /(£) is the frequency of clicks of Dl in measurements l + 1,..., k, (see (1.6 I), 
and p(L\v) has been defined right above (1.51 . In (1.161 , e = (e m )^ =1 is a suitably chosen sequence 
of positive numbers converging to 0, with sjm£ m —> 00 , as to —> 00 . An appropriate choice of e 
depends on details of the model used to describe the system S; (see Sects. and the appendix 
for simple examples). We set 


N 


"(i.fc)/ x = 


fe) = u 




(§.)■ 


(1.17) 


The following claim is obvious. 


(A) If k — l is so large that e k -i < k/ 2, where k = min |p(L|zq) — p(L|i/ 2 )|, see (1.101, then the 

"17^2 - 

sets 3^’ fc ^(e), ...,E^ k \£) are all disjoint from one another. 

Let , CfcI T !j\ -,■■■■, r )j „,), (jm < 0) denote the conditional probability of measurement 

outcomes given that the outcomes in measurements j 1 < j 2 < ... < j m < l were 

Tij 1 ,-.-,r]j rn , respectively. ( As an aside, we note that 


Mu, (£z + l I •••! \Vjl ) •••! Vjrn ) 




(1.18) 


with E(...) as in (1.111). Under quite general hypotheses on the model used to describe S, we will 
show: 

(B) The conditional probabilities of the complement of the set E^ ,fe )(£) satisfy a bound 

Mu,((S (i,fc) (£))%,,...,%•„) = l-/x w (3 ( *’' c) (e)|? ?jl ,...,J7 : , m ) < 4-z, (1-19) 

for a sequence 5 = (5 n )^=i converging to zero, as n —> 00 , uniformly in ijj 1 ,..., rjj m , for each 
fixed l. 

The bound (1.191 means that the frequencies, f^’ k \^), of long, but finite measurement pro¬ 
tocols £; + i, are almost always very close to one of the numbers p(L 11), ...,p(L\N), indepen¬ 

dently of what has been measured in earlier measurements. 

Statements (A) and (B) are analogous to the statement that, in the statistical mechanics of 
classical lattice gases, the canonical and grand canonical ensembles are equivalent in the ther¬ 
modynamic limit, independently of boundary conditions, and that the limiting Gibbs states are 
“mutually singular”, as some thermodynamic parameters, such as the chemical potential of parti¬ 
cles in a lattice gas - the analogue of v £ {1,..., N} - is varied. The mathematics used to prove 
these statements is similar to the one used to prove (B) and is borrowed from methods used in 
“hypothesis testing”, (see 6 . |20]h 
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We say that a measurable function h is an “observable at oo” (or a “pointer observable”) for 
S iff h is well defined and bounded except on sets of /z w -measure 0, and the value, h(f), of h at 
an arbitrary point £ £ S is independent of the first m components, of £, for any m < oo, 

almost surely with respect to the measure /z w , where u> is an arbitrary state of S. An example of 
an observable at oo is the limiting frequency f$°°\£) defined in (1.8). Observables at infinity span 
a commutative algebra of functions on 5, denoted by Goo [S']. Let us assume, for simplicity, that 
G^S] is finite-dimensional. Then it is the algebra of functions on a finite set 




( 1 . 20 ) 


for some finite integer N x . (More generally, Goo [5] is the algebra of bounded measurable functions 
on a compact Hausdorff space, S^; see e.g. M- But ( |1.20[ ) holds in the model described above, 
with Noo = iV; see also Section 3.) It then follows that, for an arbitrary state w of S, the measure 
Hu on S associated to u has the unique decomposition 


Mw(') — 




Pu{v)Hu(-W) 


( 1 . 21 ) 


where the measures Hui'Wi) and Hu{ m \ v 2) are mutually singular, for distinct values of V\ and i/ 2 , 
and P u {v) > 0, for all v, with Pu>{v) = 1. 

The point of these remarks is that, even without knowing, a priori, which physical quantity, 
A, of S can be measured indirectly, in a non-demolition measurement, with the help of very many 
successive direct measurements of the quantity represented by X , the commutative algebra Goo[S] 
of functions measurable at infinity tells us something about any such quantity: The operator 
A representing it corresponds to a continuous function on E^. (In the example of the system 
discussed above, A = JY, where ^ is the electron number operator associated with P.) A 
point v £ -oo is henceforth called a fact (concerning the system S). It can be reconstructed 
from very long measurement protocols. Results on “hypothesis testing” are useful to understand 
how well functions measurable at infinity can be approximated by sequences of functions of finite 
measurement protocols, and how rapidly the latter converge to the former, as the length of the 
measurement protocols tends to oo. (A detailed study of these matters will appear in forthcoming 
work; but see Sects. 2 and 3.) 

We now return to analyzing the system S introduced above; see Figure 1 and Eqs. (1.1) - (1.5). 
We consider the situation where the operator ^Y(t) slowly evolves in time t. We consider subsets, 

E ^.„(/:;£), with the property that the frequencies, f^ k fe,lfc ^(^) ) of clicks of the detector Dl in 

measurements ik — k + 1, ...,ik is within from the Born probability p(L\vf) introduced above 
(1.5), for all i = 1 If the number operator JY commutes with time evolution then the 


/z^-measure of the complement of the set 


o{k;e) ■— (kj e) 


( 1 . 22 ) 


is tiny (< pSk, where (5^)^! is as in (1.19)), for an arbitrary state co of S. If JY(t) evolves in 
time t very slowly, e.g. if 

max || ^Y(t + s) — ^Y(t) || <C 1, 

0<S<T 

then the /z^-measure of the complement of 5 p (fc;£) remains very small: 

p(k;e) c ) = 1 - /z u (S p (fc;e)) < 2 p5 k < 1, (1.23) 

provided k is large, but p is not too large. This fact enables us to introduce a measure, on 
finite “quantum trajectories” (or “histories of facts”), {v(ti) =: Vi \U = ikT} T -_ :l , with p such that 
2 pS k < 1 : 

&u,{vi ,-,Vp) := Hu(^i,...^ P {k-,e)). (1.24) 
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By (1.231, 


&w(vi, —iVp) > 1 - 2p6 k « 1. 


(1.25) 



It is interesting to characterize the measures more precisely in certain limiting regimes. 

Section 2 contains an outline of general concepts and new results concerning indirect and, 
in particular, non-demolition measurements. In Section 3, we introduce a family of concrete 
models describing systems S of the type discussed above; (see |3j [2, 4j). For these models, the 
general picture sketched above can be translated into rigorous mathematics. In Section 4, we 
study indirect measurements of quantities evolving slowly in time, and we analyze trajectories of 
“quantum jumps” (described, in a suitable limiting regime, by Markov jump processes). 

To conclude this introduction, we remark that, in many quantum systems tracked with the 
help of direct observations of probes, these probes are coherent pulses of field quanta of a wave 
medium, such as the quantized electromagnetic field, rather than electrons in a conducting chan¬ 
nel. An example is the indirect observation of the trajectory of a charged particle (P) by means of 
light scattering (E). A mathematically precise analysis of such systems is, however, quite involved 
and is postponed to a future publication. 

Acknowledgements. We thank Ph. Blanchard for numerous helpful discussions and encour¬ 
agement, and M. Bauer and D. Bernard for some correspondence at an early stage of this work 
and some discussions. M. B. thanks Jose Luis Angel Perez Garmendia for helpful discussions. J.F. 
thanks the Institut des Hautes Etudes Scientifiques for hospitality during final stages of work on 
this paper. The research of B.S. is supported in parts by the region Lorraine. 

2 Indirect retrieval of information on quantum systems 
through sequences of projective measurements 

In this section we explain in general terms how interesting information on a quantum system S 
can be retrieved from long sequences of projective measurements of a possibly very uninterest¬ 
ing looking property of S. Our analysis can be viewed as a contribution to a general theory of 
“non-demolition measurements/observations”. In this section, it is carried out within the stan¬ 
dard Hilbert-space formalism of quantum mechanics; (although using a more abstract algebraic 
formulation of quantum mechanics involving C *-algebras of operators would actually be more 
natural - see, e.g., [H]). We consider an isolated quantum system S characterized by the data 


(J^,{C/(t,s)} tiSeM , & s ), where 


• J#s is the Hilbert space of pure state vectors of S. 

• {U(t, s)}f, S £R is a family of unitary operators on M’s representing time evolution. They 
satisfy the composition rule 


U(t,s) = U(t,r)U(r,s), Vr,s,(eR. 


( 2 . 1 ) 


• ffg C 3§{Ms) is a set of self-adjoint operators on Ms representing physical quantities of S 
that can be measured/observed in projective measurements. This set has the property that 
if an operator X belongs to ffs and / is a real-valued, bounded, continuous function on R 
then f(X) belongs to Cs, t°°. 

In many examples of concrete physical systems S, the set @s is generated by only a few (often 
finitely many) operators. To simplify matters, we will assume that &s is abelian, i.e. that it is 
generated by commuting, self-adjoint, bounded operators. In this case, we may identify &s with 
the commutative algebra generated by all operators in In this paper we will usually assume 
that the spectrum, as, of the commutative algebra @s is a countable set of points. We then denote 



by 7 Tj (the orthogonal projection on iffs corresponding to) the characteristic function of the set 
{£}, where £ is an arbitrary point in the spectrum as of 0s- The operators (7r^)j gCTS satisfy 

E ^ = 7T£7T£/ = S u > 7T f . (2.2) 

Jen's 

All operators in i &s can be written as linear combinations of the projections (7r^)j gcrs . 

2.1 Calculus of frequencies and observables at infinity 

We suppose that successive projective measurements of quantities in 0s are carried out at times 
U, with ti < ti+ 1 , i = 1,2,3,.... Given a state uj(-) = Tr(p^(-)) on where p u is a 

density matrix specified at some time to, the LSW formula (for Liiders, Schwinger, Wigner; see 
[22, [22, EH, and m ) yields the a-priori probability, /r w (£i ,...,£„), for observing a sequence 
of measurement results, G £ erg, observed at times ti, i = 1,2, ...,n, corresponding to “events” 
(7r £l (ti), ...,7T£ n (i„)), where 

7Tj(i) := U{to,t)^U{t,t 0 ). (2.3) 

Indeed, the LSW formula tells us that 

M£i,-,£k) = Tr(7r ?fe (4) • • -^(iiVu^ji^i) • • -7^(4)). (2.4) 

Remark. It is well known that the LSW formula should only be used if a suitable form of “de¬ 
coherence” holds for histories of repeated projective measurements of quantities represented by 
operators in 0$', see, e.g., [15] . [14] . (where a theory of projective measurements has been devel¬ 
oped), and references given there. We will not discuss these matters here, because they do not 
play an important role in the following. 

It follows from \2A\ and from \2.2\ that 

Y /M£i,-,£fc) = for arbitrary k < oo (2.5) 

ffeGo-s 


and 

E MO = !. (2-6) 

and hence the Kolmogorov consistency criterion implies that the functional fi u extends to a prob¬ 
ability measure on the measure space (3, E), where 3 = <7g^ is an infinite Cartesian product 

of copies of the spectrum as of 0s and E is the cr-algebra generated by the cylinder sets in 3. 
The cr-algebra E is the Borel cr-algebra for the compact topological space 3. Let L°°(3, be the 
usual space of equivalence classes of bounded measurable functions with respect to the measure 
Hu- This space contains a closed subspace, denoted by 0 oo [u], of bounded measurable functions, 
h, with the property that 

h(£) = /i(£ 7 ), almost surely with respect to 

whenever G = £', except for finitely many i. Of course, L°°(E,p u ) and 0 X [w] are algebras; 0 os [u}\ 
is called the “algebra of observables at infinity” (or “algebra of pointer observables”). 

Given a function / £ L°°(3,/z w ), we define a bounded linear functional on the algebra of 
observables at infinity by setting 


m/ \h) ■■= / /modi* u (o, 


(2.7) 


for any h £ 0 oo [ui\. If / > 0, then E w (/|-) is actually a bounded positive linear functional on the 
algebra 0^ [w] and hence corresponds to a regular Borel measure p u {f |-) on the spectrum, denoted 


by Soo, of 0oo[uj\. By (2.71, E u (/|xa) = 0 if E w (/ = 1|xa) = 0, where A is an arbitrary Borel 
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subset of Soq. Thus, the measure /x u (/|-) is absolutely continuous with respect to /^(1|-) =: P u . 
Consequently, the Radon-Nikodym derivative of /i w (/(•) with respect to is well-defined, 

and we have that 




( 2 . 8 ) 


for all positive functions / £ L°°(s, /r w ). This formula can be extended to arbitrary real functions, 
/, in L°°(E:, /i w ) by writing / = /+-/_, with /+, /_ > 0, and applying p~8j ) separately to /+ and 
/_ and, subsequently to arbitrary complex-valued functions in L°°(2,Following the proof of 
[251 Th. 5.7], one can find a set TV C 2^ of measure zero and select for each / a representative 
of the Radon-Nikodym derivative such that, on the complement of N, y(^ ) is linear in /, 

positive whenever / > 0, and bounded by ||/||oo- Hence we may rewrite Eq. ( 2.8b as 


M/|i) 


/(0<W£H dP^iy), 


(2.9) 


where /x w (-|iz) is a Borel probability measure on 5, for P u - a.e. v £ E^. 

Before studying properties of the “extremal” measures 1^), v £ S^, we describe an explicit 
example of an observable at infinity: We consider a bounded measurable function / : erg 1 —> C 
and define “frequencies” /^ k \ k = 1,2, 3,..., on 5 by 

k 

/ (fe) (0 : =^E £(i)» where 

i=i 


mo ■■= 1)- (2-io) 

If the sequence of functions (/^ fe -*)^_ 1 converges n u - a.e., as k —> oo, then the limiting function, 
denoted by f^°°\ belongs to ^ooM- Sufficient conditions for convergence may be derived from 
ergodicity hypotheses. Assuming, for example, that, for arbitrary n < oo, the “decoherence 
condition” 

M£i.-,&,,e (>n) ) = M u4 (>n) ), (2-ii) 

Cl,...,Jr.ecr s 


holds, the ergodic theorem implies that converges a.e., as k —> oo, to an element in 
^ooM- For an autonomous system S' (i.e., for one with U(t,s) = U(t— s), for all t, s £ R) with 
the additional property that there is “perfect decoherence”, for arbitrary sequences of projective 
measurements of quantities represented by operators in Eq. (2.111 holds. 


2.2 Hypothesis testing and emergence of facts 

Let / : (Tj 1 -t C be a function depending on finitely many measurement outcomes, and let /^ 
denote the frequency associated with /, as defined in (2.101. We decompose Soo into the level sets 
of the limiting function /(°°): 

N 

Soo = IJ ( 2 - 12 ) 

a=l 


where is a level set of /(°°). We henceforth assume that 

min ai ^ Q2 |/ (oo) | H «i - / (oo) | h «2 | =: k > 0. 

We define 

'■= J d/j, u (-\u)dP u (v). 


(2.13) 


(2.14) 


Note that the measures dfp^ are mutually singular for different values of a £ {1,...,IV}. We 
introduce the “fluctuation variables” 


<^ fc) := Vk (/ (fc) - m^j , where m a := / (oo) (0Ie 


(2.15) 
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Our analysis is based on the assumption that (d/z“, </>&) satisfies a suitable Central Limit Theorem 
mmm, for all a = 1,2 We consider the generating functions, Fi k \h), of connected 
correlations of fluctuation variables defined by 


Fi k \h) :=\n(JjvG{Z)e h *™®y 


(2.16) 


where h £ R. We assume that, for every a = 1, there is an open interval (JcK containing 

the origin such that the sequence (Fi k \h)) k T 1 of real-valued functions, together with their first, 
second and third derivatives in h converge to bounded, continuous functions of h £ U, as k —> oo. 
Furthermore, we assume that - to mention a concrete example of a suitable hypothesis - 

• ^(OMO, 

• -S' 7a > 0, and 

• £?Fi k) (h) -»■ 0, for all h £ U, 

(k) 

as k — > oo, for all a = 1,2, ...,7V. Then the fluctuation variables (jy a ; approach Gaussian random 
variables, with variance "f a centered at 0. This is the contents of the Central Limit Theorem. 
The hypotheses formulated above can be derived from suitable cluster properties of the cumulants 


^h,k(fji> fj m )i for m = 3 and all h G U, 

where k denotes an expectation with respect to the probability measure 

exp[-Fl fc )(h)] dfiH)- 


To be somewhat more precise, these cumulants should be assumed to be summable in ... ,j m , 
for every fixed j -\, with bounds that are uniform in k. 

We introduce subsets 

4 fe) fe):={|l \f (k \i)~m a \<e k }, (2.17) 

see Eq. (2.151, where is a sequence of positive numbers converging to 0, with e k ■ Vk — > oo, 

as k — > oo. 

The following theorem is an immediate consequence of (2.131 and of the Central Limit Theorem. 
Theorem 2.2.1. 1. If k is so large that e k < § then 


s^GOn h«(s) = i 


(2.18) 


for ai^a 2 - 

2. The measure of the complement of the set U Q si^(£) tends to 0, as k tends to oo; i.e., 

„ u ((U a 5W(£))j^0, (2.19) 

(k) 

as k oo, at a speed determined by the speed of convergence of the fluctuation variables < jy a to 
Gaussian random variables. 

This theorem tells us that the value, a, of /(°°) can be inferred from long, but finite measure¬ 
ment protocols, with an error that tends to 0, as k approaches oo. The value of /(°°) will 
henceforth be called a “fact”. Apparently, long measurement protocols reveal facts about the state 
of the system, and the chance of being mistaken about which fact is emerging tends to 0, as the 
length k of the measurement protocol tends to oo. 
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The a-priori probability of observing an arbitrary measurement protocol ( < ' k> corresponding to 
the value a of /(°°) approaches the value 


:= J dP u (y) 


( 2 . 20 ) 


This is Born’s Rule. 

Remark. With the frequency / ^ we can associate an operator, A^ k \ acting on the space of 
density matrices and defined by 

1 k 

A(k \p ) : = £ /(&>-jO+j-i) nj+i-Atj+i-i) ■ ■ ■ nAtj) p 

i =1 £j,...,£j+i-i 


where p is an arbitrary density matrix. Viewed as an operator acting on density matrices, an 
“observable at infinity” is obtained by letting k tend to oo in the above expression. In concrete 
examples of quantum systems S, it is often possible to identify an “observable at infinity” with 
an explicit operator, jY , (with \<j(<jY)\ = TV), acting on the Hilbert space M’s of the system 
that has the property that it commutes with time evolution; i.e., M(t) = U(to,t)MU(t,to) = 
M independent of time t. One then speaks of a non-demolition measurement of the quantity 
represented by the operator M. 

To conclude this section, we remark that, in practice, the algebra of observables at infinity 
is most often trivial (i.e., it consists of multiples of the identity), because most “facts” about a 
concrete quantum system S tend to actually vary with time, and the revelation of such facts, 
while usually accomplished through long sequences of indirect projective measurements, does not 
correspond to a “non-demolition measurement”, in the strict sense of the word. In this situation, 
it remains meaningful to attempt to construct a suitable decomposition of the support of the 
measure dp u into subsets corresponding to “facts” that are only valid (with very high probability) 
during some long, but finite interval of time. To be concrete, we define functions, j( fc ’ fc + r ') ) where 
r = 0,1, 2,... is a coarse-grained “time variable”, as follows: 


k-\-r 


f (k ' k+r) :=l £ 


j-k+l 


see (2.101. We introduce sets 

~i k ’ k+r \e) := U | |/(W)(|) -m a | < e r }, 


( 2 . 21 ) 


( 2 . 22 ) 


where m a is chosen suitably - it replaces the quantity |s~, ol = 1,-..,TV in Eq. (2.17l 

and the sequence (e n )^ =1 tends to 0, as n —> oo, as above. We assume that 


Clearly, 

if r is so large that e r < |. 


min ai ^ a2 | m ai - m a2 \ > k > 0. 

E^’ k+r \e) nE^’ k+r \e) = 9 

Next, we introduce the sets 

3a,. , a Jr;e) := f] E^~ r ^(e). 

j=i 


We say that the sequence (a±, ...,a p ) represents a “history of (plausible) facts” iff, for a suit¬ 
able choice of positive integers r and p > 1, the set U aii ... iQ; S aii .... |Ce (r;e) has essentially full 
Pu- measure, i.e., 

|J E aii .... >ap (r;e)) > 1- <W, ( 2 -23) 

ai,...,Q!p 


12 




where S p r is so small that, for all practical purposes, it is indistinguishable from 0. If (2.23) holds 


then, with probability very close to 1, one of the histories (ai,...,a p ) will take place, and the 
probability to observe the history (oq,... ,a p ) is predicted to be given by 


(^ 1 ; • ■ ■ 5 ®-p) d’u (- 


, (r;s)), with ^(ai,...,a p ) « 1, 


(2.24) 


which is a generalization of Born’s Rule, see (2.201. 

In Section [4j we consider “interesting” physical quantities of a quantum system S that vary 
very slowly in time and can be measured indirectly through long, but strictly finite sequences of 
projective measurements of some rather “uninteresting” quantities. 

To return to the analogy between the quantum systems discussed above and one-dimensional 
lattice gases (or spin chains) alluded to in the introduction, we point out that “time” in the 
quantum systems becomes “space” in the lattice gases , measurement protocols, £, are analogous to 
configurations of (various species of) particles, the measures p u are analogous to Gibbs equilibrium 
measures of lattice gases, the frequencies /(°°) are analogous to particle densities, with a playing 
the role of a chemical potential, Theorem |2.2.1| is analogous to results concerning the equivalence 
of the canonical and the grand-canonical ensemble in the thermodynamic limit, k —> oo, of lattice 
gases. In the theory of lattice gases, it is well known that the Gibbs measures corresponding 
to different particle densities are mutually singular, (i.e., have disjoint supports). The material 
discussed between Eqs. (2.211 and (2.241 would correspond to a lattice gas with a space-dependent 
chemical potential and hence a spatially variable particle density. 


2.3 Exchangeable measures on the space of measurement protocols 

We conclude Section 2 by considering a special case of the general theory of non-demolition 
measurements developed in the last two subsections. We assume that the probabilities of results 
of successive projective measurements of observables in are independent of the order in which 
these results are observed. To be more concrete, we imagine that, at each time of measurement, 
the quantities corresponding to the operators in 0$ are measured with the help of some “probe” 
temporarily interacting with the system, and that the probes used at different measurement times 
are different from each other and are prepared and measured independently of each other. A 
specific model of this type has been described in the Introduction and will be studied in more 
detail in the next section. Mathematically, the above assumption is translated into the following 
property of the measures p u \ 


/4j(£tt(i)) ■■■)^ 7 r(fc)) is independent of the permutation 7r, 


(2.25) 


for arbitrary permutations, 7r, of for all k = 2,3,... . Measures satisfying (2.251 are 

called “exchangeable”. Such measures have been characterized completely by de Finetti and his 
followers laum who have shown that such measures can be written as convex combinations of 
product measures. Thus, if is exchangeable then there exists a measure space Soo, a measure 
dPu on Sqq and probability distributions p(ji') on a$, v G S^, such that 


P k 

M£i> ••■>&)= / dP u (y) 

Js °° j=i 


(2.26) 


This equation is a concrete instance of the general representation (2.91. It follows from Kol¬ 
mogorov’s zero-one law that Soo appearing in the decomposition is indeed the spectrum of the tail 
algebra. For exchangeable measures, we introduce the frequencies 


■= and := -P(£W))- 


1 


(* 0 / 


c (G / 


o =i 
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Using the Law of Large Numbers and the Central Limit Theorem in the form due to Berry and 
Esseen, see, e.g., it is straightforward to prove an explicit (quantitative) version of Theorem 

2 . 2 . 1 . 

In subsequent sections, we study families of models of systems S giving rise to measures 
that are exchangeable for arbitrary states ui of S. Such models have the advantage of being quite 
easy to analyze mathematically. 


3 A family of simple models of quantum systems 


3.1 Description of non-demolition measurements 

A general mathematical formalism to describe the statistics of measurement outcomes and the 
reduced evolution in an indirect measurement was developed by Kraus [21], (see also [19])• Here 
we briefly recall some elements of this formalism. We consider a system S composed of a subsystem 
P of interest to an experimentalist and a subsystem E consisting of a measurement device. We 
suppose that M’s = Mp®Ms- The reduced time evolution of the subsystem P, conditioned upon 
a measurement event p C as, is encoded in a completely positive map <h*(p) acting on Z$(Mp). 
To give rise to a probability measure on a space of measurement events this family of maps must 
be countably additive, (in particular, for disjoint sets pi and p 2 , <h*(piUp 2 ) = (pi) + 4>*(p 2 )), 
and <&*(o\s) must be the identity map on 33{Mp), (i.e., $*(cts)[ 11] = 1). The maps $*(•) describe 
the evolution of “observables” of the subsystem P in the presence of measurement events (i.e., 
events happening in the device E)\ the evolution of states of P is described by dual maps, $(p) = 

For simplicity we assume that as is finite and that dim(Mp) < oo. Hence $*(•) is completely 
determined by the maps 4>*({£}) = : 4>*£, with £ £ as', in particular, <F*(crs) = ^*{- For 

a given state, w(-) = Tr(p(°) (•)), of P, the probability to observe a measurement result £ is 
given by w(< 1>*£[11]), and the reduced non-normalized posterior state of P, given that £ has been 
observed, is For a sequence of independent identical probes, as discussed in Section 

the probability to observe the measurement protocol = (£ 1; ...,£*.) is given by 


2.3 


/^(£ (fe) ) = w($* Sl O ■ • • o $*£*[!]), 


(3.1) 


see (2.4). The corresponding reduced state, p( fe )(£U)), of P is then given by 

O <l>c, >H'>] 


,(k) 


(£«) = 






Tr(% o---o$ ?1 [p(0)])' 


(3.2) 


Our goal is to find the tail algebras for all such i.i.d measurement processes, i.e., to determine 
the “facts” emerging, in a given experiment, from projective observations of very (infinitely) many 
identical, independent probes that have interacted with the subsystem P. In this paper we wish 
to determine these tail algebras in the setting described in Sections 1 and 2. 

In a non-demolition measurement, the mere presence of a measurement device does not affect 
the probabilities of sequences of measurement results. Mathematically, this translates into the 
condition 

'y ' Puijii, ■ ■ •, £,j— i) £.7+1) • • • > £fc) = • ■ ■, i, £j+i) • • • i £fc)- (3-3) 

fj6crs 


Upon inspection of Eq. (3.11, we see that a suitable assumption that guarantees that this condition 
holds is 

<I>*£ o <!>*£/ = o for all £, £' e a s . (3.4) 

(A more detailed discussion of such hypotheses will be presented elsewhere). Condition ( |3.4| ) 
immediately implies that the measures are exchangeable, see (2.251, and hence has a de 
Finetti decomposition; see (2.26). Assuming that condition (3.41 is satisfied, we can fully describe 
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the “algebra of pointer observables” under the technical assumption that $(< 75 ) has a faithful 
stationary state. In this particular case, we show in Lemma |5. 3. 2| that 


[$* £ [l],$* e [ll]] = 0, for all £, £'£ <7 S . 


(3.5) 


We denote by {11^|II^ £ v £ Sqo}, the joint spectral projections of the commuting family 

{$*£[l]}£eo-s- The conditional probabilities p(£\v) appearing in the decomposition (2.261 are the 
eigenvalues of <I>*£ corresponding to the eigenprojections n^, i.e., 


$*#] = E 




and the probabilities p(*|z/) on as are mutually distinct. Eq. (2.26) then takes the form 


(3.6) 


MO = E w ( n 0 M|H with MM-MIO = YlpfaW)- 

i/eHoo j =1 


(3.7) 


The tail algebra is isomorphic to the commutative algebra generated by the orthogonal projections 
i.e., by all operators of the form Yh v /(i / )n i/ , where / : S,*, —> R is an arbitrary 
bounded measurable function on S,*,. In p!]. Maassen and Kiimmerer have considered the special 
case where 

= C ( pC* 0 (3.8) 

Cj being the operators on that appear in a Kraus decomposition, 


*(°s)\p] = E mm E c l c i = a - 

{Go’s { 60 's 


(3.9) 


A special feature of this choice is that each map <f>£ maps pure states to pure states. The family 
of operators satisfies (3.4 ( and (3.6 1 , provided that C^ C£(i/)n„, with \c.(v)\ yJ 

|c.(z/)| if v z/. A straightforward calculation then shows that 

$*j[l] = E MHM where p(£\u) = |c 5 (^)| 2 . 


The normalization condition in (3.91 ensures that p(-\u) is a probability distribution. In a non¬ 


demolition experiment of the kind described in the Introduction, the coefficients c^{y) can be 
interpreted as transition amplitudes for a probe (always prepared in the same initial state, ipo) on 
which a projective measurement is performed after it has interacted with P, with the measurement 
result corresponding to a final state ip£. The Hamiltonian describing the time evolution of the 
probe interacting with the system P depends on the point v £ i.e., the joint Hamiltonian of 

the system P and a single probe has the form H = n„ ® H v , which leads to a joint time 

evolution (over one measurement cycle) of the form 


u = Y n „®U V . 


(3.10) 


i/GE 


The coefficients c^{y) are then given by c^{v) = {ip^,U u ip 0 ), where ip 0 is the initial state of the 
probe and ip£ is its state after a projective measurement on the probe sub-system has yielded the 
result £. 


3.2 Emergence of facts in non-demolition measurements 


In this section we proceed with the analysis of the measures given in Eqs. (2.261 and (3.71, 
under the additional assumption that the spectrum Soo consists of finitely many points and, hence, 
can be identified with a subset of R. We assume that as is countable and allow 3@p to be infinite 
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dimensional. It is well known that the product measures appearing in Eqs. (2.261 and (3.71 
satisfy large deviation estimates. For every fixed £ £ < 75 , we introduce the random variable 


J?E) = jf fe) (£) := \#{j e {!,••• ,fc} | £i = £} 


(3.11) 


on (2, E). Given v £ 2^, the empirical measures satisfy large deviation estimates 

below. Before stating this 
iminary facts. 


3.2.1 


and converge to X){eo- s P(£\ l> )^£ exponentially fast; see Theorem 
result more precisely, we introduce some notation and recall some pre 

We denote by ^ft{a s) the polish space of probability measures on a s', (we note, in passing, that 
as is a polish space, because every countable compact Hausdorff space is polish). Convergence 
in „-#( as) is equivalent to weak convergence, i.e., to convergence against bounded continuous 
functions on as- Let p = p{£,)^ € ■3C(as) be an arbitrary probability measure on as- We 

denote by 

Vv = E p(£ W) S Z (3-1.2) 

Jetr s 

the probability measure on as with density p(£ \v). The relative entropy of the probability measure 
p with respect to the measure p u is defined as 


J pAp) ■■= 



^) ln ®> 


ifp < Pu, 
otherwise. 


(3.13) 


With an arbitrary subset I\ C {as) of probability measures on as we associate the quantity 


J pA k ) ■= inf I p Ap)- 

p€K 


(3.14) 


It is easy to see that I Pv is a convex function of p and that I Pv {p) > 0, for every p £ ^#( as ), with 
equality iff p = p v . The following theorem is well known; (see, e.g., El ED])- 

Theorem 3.2.1 (Boltzmann, Sanov). For any closed subset K of ^/K{as) not containing p v , one 
has that I Pu ( K) > 0 , and there is a constant Ck > 0 such that 


I E 4' fc) 0£)% G K}\u) < C K exp(— (K)/2) 


€'eo-s 


(3.15) 


for all k £ N. 


The construction leading to Eq. (3.71 ensures that the probability distributions p u and p u 


are distinct, and hence the measures /x w (- \v) and p- u {-\v') are mutually singular. We can then 
decompose S into a collection (2^) l/e H oo of measurable subsets of 2 such that fj, u {S v i\i>) = 8 v y. 
The proof of the mutual singularity of the product measures /r w (-|z/) and p w {- \v') is straightforward 
(involving a lim inf of sequences of sets). In order to keep this paper as self-contained as possible, 
we sketch it in an appendix, Sect. 5. As already discussed in the Introduction and in Section 
2, Theorem 3.2.1 enables an experimentalist in the lab to measure properties of the system P 
corresponding to functions on 2 ^ indirectly, and hence to determine v £ 2 m , by recording 
empirical frequencies of long sequences of outcomes of projective measurements of (possibly very 
uninteresting) quantities in (3s- 


3.3 Decoherence and “purification” 

In the previous subsection, we have shown that, in non-demolition experiments, facts v £ 2^ 
(identified with a subset of R) emerge exponentially fast. In this subsection, we are interested in 
studying the time evolution of the reduced density matrix describing the state of P, as projective 
measurements of quantities in @s are carried out, with results £i,...,£fc belonging to as, with 
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k = 1,2,3,.... In particular, we will prove that, under suitable hypotheses, the reduced density 
matrix of P has a limit, as fc —> +oo. 

At the level of generality chosen at the beginning of Section |3.1| the limit of these reduced 


density matrices may not exist, because (3.6 1 does not sufficiently constrain the time evolution 
restricted to the range of !!„. We avoid this difficulty by assuming that the time evolution of the 


system consisting of P and one probe has the form given in Eq. (3.101. Thus, in this section we 


consider maps <f>e given by Eq. (3.8 1 , with (cf. text below Eq. (3.91) 


C t= 


In summing over we tacitly assume that the measures p v (see (3.12)) are mutually distinct. 
The dynamics, constrained in this way, has the property that the maps 4>*£ act on operators of 
the form n„(-)n„/ as multiplication operators, and we have that 


n ^( $ «[p]) n ^' = c^{v)Ii. v pU. v >c^{u') 


(3.16) 


for an arbitrary density matrix p on dFFp. Next, we state and prove a result of |24| rediscovered in 
[3]. Our proof is significantly different from those in (241 [3j; (it shows that the result is a nearly 
immediate corollary of Theorem 3.2.11. In the following, || • || denotes either the operator norm or 


the trace-norm, the statements being correct for both choices. 

Theorem 3.3.1. For v ^ v', 

||n^ fe )(^)n^||-to, Hv-a.8., (3.17) 

as k —* oo, where p( fe )(£( fc )) has been defined in (3.2). Furthermore, there exists a random variable 


0 : 


such that 


,(*0 


(£«) _ 


iW 0 ) n e 


p u - a.S; 


Tr(n e p(°)lle)) 

as k —> oo. The probability that 0 = v is equal to P u {y) = Tr(p^II„). 


Proof. Eq. (3.16l yields a recurrence relation for the density matrix p^ k \ 

Jk-i) ( Ak-i) 


~ " E |c ? »| 2 Tr(n^(fc-D^-D))’ 




(3.18) 


(3.19) 


which, by iteration, implies that 


iW fc )(!«)iv 


i 

Pui{£, i, £&) 


k 

n y/ 9 (0) n y / n c «» c «>')- 

i =1 


(3.20) 


Multiplying both sides by by ...,£fc), taking norms, and summing over all possible £i, 

we find that 


E||n^( fc )(e (fc) )n^|| < ||iW°)iv||(^ ]c|(i/)| \c^)\) k < \\n„pWn„,\\(s„„,) k , (3.21) 

JetT S 


for some 0 < <5< 1. The factor (<5„„') fc on the right side is obtained from the Cauchy-Schwarz 
inequality and the non-degeneracy assumption, p v ^ p u * for v ^ v'. The exponential decay of the 
expected value in Eq. (3.211 implies almost sure convergence - a consequence of the Borel-Cantelli 
lemma. 


17 



























To prove (3.181, we observe that from the fact that the measures //„(• \v) are mutually singular, 
as v ranges over one may deduce that the sets (introduced at the end of Section 3.21 have 
the property that, for v A i/, 


(e (fc) )n„ I|xs„, ->0, as 00 , 


(3.22) 


Hu- a - e -- This property is proven in the appendix (see subsection 5.1.2), where we will also show 
that 

||TV(n^°))n^ fc )(^ fe) )n,-n^ 0 )n I/ ||xH„( 0 ^o, ^-a.e., (3.23) 

as k —» 00 . Let 0 be the random variable defined by 

©(C) = 55 *=„(£> v. 


Using (3.21), (3.22), (3.23), the triangle inequality and the fact that /i u (U„s„) = 1, one shows 
that 


4*0 


(C W ) - 


ne(j)P (o) n 0(?) 


Tr(n e(? )p(°)ne( 0 )) 


converges to zero, fiui- a.e.. 


□ 


4 Indirect measurements of physical quantities varying slowly 
in time 

4.1 A perturbative approach to indirect measurements 

At the end of Section |2.2[ we outlined some general ideas concerning indirect measurements or 


observations of physical quantities that evolve slowly in time; see Eqs. (2.211 through (2.24). In 


this section, we intend to add some mathematical precision to those ideas by analyzing indirect 
measurements of a physical quantity evolving slowly in time, for a class of simple models. For this 
purpose, we develop a perturbative approach to the theory of indirect measurements. We consider 
models whose time-evolutions are perturbatively close to the ones considered in our analysis of 
non-demolition measurements presented in Section |3.1| As in Section |3.1| we consider models of a 
quantum system S that is the composition of a subsystem P of primary interest with a subsystem 
E consisting of equipment used to observe P. We are interested in describing observations of 
physical properties of the subsystem P, also called “facts”, using the experimental equipment 
described by E. For this purpose, successive projective measurements of some quantities referring 
to E are carried out. They yield a sequence of reduced states, p ( ' k \A' k ' > )i on the algebra, 
of bounded operators on depending on measurement protocols k} k> = (£i, G {<Js) xk 

of arbitrary length k < oo, (with £j the values of a quantity referring to E measured in the j th 
projective measurement). 

The reduced state of P, after k projective measurements carried out on E, is assumed to be 
given, recursively, by the density matrix 




$fV fe 1) (£ (fc 1) )] 


Tr ($f ) [p( fe - 1 )(^ (fc - 1) )]) : 


(4.1) 


where, for all £, the evolution map <1>^ ' ) only depends on the first k measurement results = 

(£i, ...,£&), for all k £ N. As before, we assume that each map is the dual of a completely 
positive map and that 

55 Tr(4>f ) [p])=Tr(p) (4.2) 

Sfeeo-s 
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for any trace-class operator p £ i^(j#pj. The probability of the measurement protocol f} k \ 
assuming that the system has been prepared in the state u> corresponding to a density matrix p^°\ 
is given by 


fc (ei,,.,a) = h(fo...o ^V 0 ) d. 


(4.3) 


In contrast to the assumptions required in Section |3.1[ we allow the commutator appearing 


in Eq. (3.51 to be non-zero, but small, (and the system is usually not autonomous).We compare 
the evolution described in (4.11 to one considered in our analysis of non-demolition measurements 
which it is assumed to be close to. We add a tilde to all quantities referring to non-demolition 
measurements, as treated in Section [3.3| namely, we will write 


$t[p] : = ^(v)n v pH vl c^(v'), 

/M£l,-,£n) := Tr($^„ o...o<%[p(°)]). 


(4.4) 

(4.5) 


To avoid confusions, the trace norm of a trace-class operator p is denoted by ||p||, and the operator 
norm of an operator A in is denoted by ||A|| op . The operator norm of a map 4^ (= <f>t 

or 4>^), viewed as an operator on the Banach space of trace-class operators, is denoted by ||4/||. 
In this section, we always assume that erg and Soo are finite point sets. We propose to study the 
dynamics of the reduced system P under the following assumptions: 

Assumption 4.1.1. We assume that there exist constants d\ £ [0,1) and £ (di,l] such that, 
for all n £ N, for every £ £ as and every f with 

(i) ||$f -^11 <di||$£||, 

(ii) Tr(<E <^p) > cfell^dl* for all density matrices p on Jftp. 

The first assumption concerns the smallness of the difference between the actual evolution and one 
corresponding to a non-demolition measurement, while the second assumption implies that the 
a-priori probability of any outcome of a projective measurement of a quantity referring to E (used 
in a non-demolition measurement of a quantity of P) never vanishes. Items (i) and (ii) imply that 

T^($f [p])>(d i -di)[|%||, (4.6) 

= :d 


for an arbitrary density matrix p and any k. Without loss of generality, we also assume that 
d = d 2 — d\ < 1. Before stating the main results of this section, we give some examples of maps 
satisfying Assumption |4.1.1[ 

Example 4.1.2 (illustrating (i)). Let (P^)fcLi he a sequence of Hamiltonians with ||P( fc )|| op < 
d\/2, for all k. Then (i) holds for a perturbation dynamics = <l>e fc exp(— iAd H (k )). 

Example 4.1.3 (illustratig (ii)). Let ffy = Vc id + T^, where v j is a probability distribution on 
ag and is a family of maps satisfying Tj = 0. Then 4>t satisfies (ii) with the constant 


^2 


sup 

Jeers 


1 - 2 


IlTdl \ 

«c + irw 


provided > ||T|||, uniformly in If the ratio of ||T^|| and v £ is small, for all then is close 
to 1. 

In studying the first example we use the inequality |1 — e zx \ < |a:| and conclude that 

ll$f° -<%J| < ||^eJ|||l-exp(-*Ad ffffc ))|| 

<2|[^J||P«[| op , 
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where we have used that ||Ad#|| < 2||i4|| op . To treat the second example, we remark that (when 
viewed as a map acting on the space of trace-class operators) satisfies the inequality 

(^-||T c ||)-td<^<(^ + ||Tc||)-id. 


Hence we have that 


Tr(<l 5 p) >V£- || || : 

> 11 ^ 11 ( 1-2 


^ + iiT ? ir 


2||Tdl 


for an arbitrary density matrix p. 


4.2 Trajectories of quantum jumps on Soo 

We denote by 

£&{■&?) ( 4 - 7 ) 

i/EHoo 

the operator representing the physical quantity of P that we wish to measure indirectly. (It is 
assumed here that Sqo is a subset of R and that v is a function on that separates points of 
In the example discussed in the Introduction, the operator JY is the number operator counting the 
number of electrons in the component P of the quantum dot P close to the conducting channel.) 
The time evolution of YY is assumed to be non-trivial under the dynamics corresponding to the 
maps <Ic fc) . Nevertheless, one may hope that the measurement protocols can be used to track 
the values of YY during fairly long, but finite intervals of time (on which the value of T is constant 
with very high probability), because the dynamics determined by the maps is assumed to 
be close to the “non-demolition dynamics” studied in Section |3.3[ see Assumption |4.1.1| For each 
k £ N, we introduce an “estimator” ^Y^ k ' k+r ^ whose value correctly predicts the outcome of 
a direct measurement of ^Y at a time « k + r with high probability, provided the outcomes, 
(£fc_l_i, ..., £fc+ r ), of r projective probe measurements are known, for some “time constant” rtl: 
For every fixed £ £ erg, we set 


/?’ fc+r) (0 := \*{j e {k + 1, ■ • • , k + r} | & = e}, 

and we define the estimator Y^( k,k+r \t;) by 

^ (fc ’ fe+r) ( 0 := argmin / p „ ( £ jf’ fc+r) (0 «5 5 ), 

Z'GS OO tr . 


(4.8) 


(4.9) 


for any fixed choice of argmin in case there does not exist a unique minimizer. We remind the 
reader that the measure p v has been defined in (3.121, and that the coefficients p(£\u) are those that 
have been introduced in the “non-demolition model” discussed in Section 3.3 i.e., p(£|i/) = |cj(^)| 2 . 
The level sets of the estimator ^Y( k < k+r '> (£), 


:= {C € S | ^ k ’ k+r \0 = u), (4.10) 

yield a partition of the space S of infinitely long measurement protocols. The probability of an 
error in the prediction by the estimator of the outcome of a direct measurement of the value of 
,/V at time k + r is then given by 

e (fc,fc+r)( J/ ) Tr((l-n„)p (fc+r) (| (fc+r) )) jU„(| (fc+r) ). (4.11) 
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The error probability for the “non-demolition model” is denoted by eS k ' k+T '(y). In the context 
of the “non-demolition model”, we are facing a classical parameter estimation problem that is 
equivalent to a problem concerning multinomial tests [T8| and hence can be solved using large 
deviation inequalities, as in Theorem |3.2.1| In particular, it is possible to characterize the speed 
of convergence to the values v on the level sets of the estimator, for the “non-demolition model”. 

The hypotheses underlying the results in this subsection are the ones summarized in Assump¬ 
tion |4.1.1| Our first result is the following lemma. 

Lemma 4.2.1 (Estimation Fidelity). There are constants a £ (0,1) and C > 0 such that, for all 
k, r £ N, the following bounds hold uniformly in p 1 ' 0 ' 1 : 


*+-)=„) -ETr(n v/ 5W)| < 2 Car, 

(4.12) 

pu(^^/p(k,k + r-) = „) ETr(n„pi i) 

<2 Ca r + , 1 d- r -\ 

d~ [ — 1 

(4.13) 

g(*,fc+0( I/ ) < Ca r , 


(4.14) 

e (M+r)(„) < Ca r + 

dl j-r-l 

(4.15) 

d-1- ^ ’ 


for all v £ Sqo. 

The proof of Lemma |4.2.1| is not particularly difficult but lengthy, and we defer it to the appendix; 
(see Section 5.21. We note that, in the “non-demolition model”, 

ETr(n,pW) = Tr(lW°)), 

and hence ETr(n„/5( fe i) can be replaced by Tr(n„//°i) in (4.121. The basic ideas of the proof of 
Lemma urn are quite straightforward. First, one controls the “non-demolition dynamics” using 

large deviation estimates that lead to (4.12) and (|4.14|. Then one controls the difference between 

' Ik)’ (k— 1 ) 

the true dynamics, as given by the composition o •• • of completely positive maps, 

and the “non-demolition dynamics”, using simple perturbative estimates. This leads to the bounds 
given in (4.13) and (4.15). The constant a turns out to be directly related to the infimum over 


v ^ v' of the relative entropies 7 p „(?v); see Section 5.2 


We deduce from the definition of the error probability e^ k,k+r ' ) {v) that 


| min Tr((l - n„)p( fe +”)(,0) > a} n i 


c (k,k+r) 


J i k + r ) =i/ j — 


H 


A 


(4.16) 


for any A > 0. Since, by definition, the union of the sets S y(k,k+r)— v is the entire space S, we 
conclude that 


Mu,({e | nun Tr((l - Y +r \t)) > a}) < £ 

“°° i/GS oo 


F (fc,/c+r 


iy) 


A 


(4.17) 


uniformly in p(°\ Combining this with Eq. 

Theorem 14.2.21 

Theorem 4.2.2 (Jump process). Let e £ (0,1]. If r is large enough and if d\ is small enough, 
then 


(4.151, we obtain the first part of our main result, 





max Tr(II^pi fc+r i(^)) > 1 — el) > 1 — e, 


(4.18) 


for all k > 0, uniformly with respect to the initial condition p^. Furthermore, 

| 3z,£~ 00 : \\p {k+r \§)-H,Y+ r \m v \\ <e}) > 1-e, 


(4.19) 


uniformly with respect to p^. 
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Theorem |4.2.2| says that if the dynamics of the system S is very close to a “non-demolition 
dynamics” then the evolution of the reduced state of P is very close to one described by a stochastic 
process on the space of density matrices commuting with the operator M, uniformly in the initial 
state In other words, the dynamics of the reduced states of P is encoded into a jump process 
on the spectrum, Sqo, of M. Unfortunately, we are unable to determine the transition rates of 
this jump process for general models. However, in certain limiting regimes, it is given by a Markov 
chain with explicit transition probabilities. Figure 2, below, refers to a simple model of this kind. 

M M 

- > 

< Ai * < Ai * 

< -> 

1 cycle of duration A 2 


In Figure 2, the duration, A 2 , of a full measurement cycle is much larger than the length, Ai, 
of the time interval during which projective measurements in E are carried out. In time intervals 
of length A 2 — Ai, the evolution of the system S is unitary; the corresponding unitary propagator 
being given by 

° 5 

for some Hamiltonian Hs acting on Ms = Mp ® Me of the form 

In time intervals of length Ai, the unitary evolution is “interrupted” by M projective measurements 
of a quantity of E represented by a self-adjoint operator X acting on M E \ (we refer to the 
Introduction for a simple, concrete model; and to m for an outline of the theory of projective 
measurements). The times of measurement of X, during the (n + l) st cycle, are approximately 
given by 

nX2 + M - l Al ’ j = 0,- 1. 

If M is sufficiently large and Ai is small enough, we can use Theorem |4.2.2| to approximately 
calculate the transition probabilities for jumps on the state space 5^ = {1,...,7V} during each 
cycle. It is not hard to show that, in the limit where first M —> 00 and then Ai —> 0 the transition 
probabilities for jumps from v to v' (with v, 1 / £ Sc») approach the ones of a Markov chain whose 
transition probabilities are given by 

Tr(n„e iX2H rIl„,e- iX2H Pl U). (4.20) 

To conclude this section, we present the proof of the second inequality in Theorem |4. 2. 2| 

Proof. (Second inequality of Theorem \4.2.2\ ) We first remark that 

EHn^+^n^H = Y, ll n ,p ( ' s+r) (! (fc+r) )n^|l M| (fe+r) ) 

= ]T ||n^ +fc) o...o$( fc + 1 )[pW(^ fe ))]n^|| M£ (fe) ) 

£(fc + r) 

< 0W) r + jpr— 
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where we have used (3.21 


4.1. 1| (see also Appendix 


and a simple perturbative estimate based on item (i) of Assumption 
5.2 for similar, more detailed calculations). The above inequality and 

(4.21) 


(4.18|) enable us to choose r so large and di so small (depending on r and c^) that 


Y, E[||n^+ r )lV||J <e 2 /4, 

v^v' 


and such that (4.18) is fulfilled, with £ replaced by e/2. Introducing the sets 
M,+r := | max Tr(n j/ p (fc+r) (£)) < 1 - e/2 j, 


and 


h+r := {{ | Y W^P (k+r \mA > e/2}, 


(4.22) 


(4.23) 


v^v' 


for all k > 0, it follows from ( |4.18 1 and (4.21) that fiuii^k+r) < e/2 and n u {3§k+r) < e/2, for all 
k > 0. Let £ ^ s^k+r U £$k+r- Eq. (4.221 shows that there exists some vq £ 5^ such that 


Y Tr(n,p( fc+r )(£)) = Y l|n,p (fc+r) IM < e/2, 


while Eq. (4.231 implies that 

Y l|n,p (fc+r) (£)n^||<e/ 2 . 

v^v' 

Gathering (4.24)-(4.25), we finally deduce that ||p( fc + r )(£) — n i/0 p( fc+r ') (£)!!„„ || < e. 


(4.24) 

(4.25) 
□ 


5 Appendix 

In this appendix we present all proofs that have been omitted in the main body of the text. 


5.1 Proofs of results in Section [3] 

We begin this subsection by explaining why the measures /x w (- \v) and \v') are mutually sin¬ 
gular, unless v = v'. Afterwards, we complete the proof of Lemma [3.3.1 


5.1.1 Comments concerning the mutual singularity of the product measures /i ul (-\i/) 

Since the measures ^ p„/ unless v = z/, we choose open sets (defined by the metric on 
..#(as )) centered on such that n = 0 unless v = v'. For every v £ 5 ^ we define 
s/ n) C Soo as the inverse image of the set under the empirical measure • It 

directly follows from ( |3.15[ ) that 

^(5^?V) < M( 3 k"' ) ) c | I/ ) < CW„ exp(-n/ p „((^) c )/2) (5.1) 

if z/ z/, and hence Borel-Cantelli’s Lemma implies that 

/x w (limsup(S[j l) ) c |j') = 0 = /K^(limsup(s[ / ? ) )|i/) = 0. 

We define S„ = }limsup(s/”- ) ) c } = liminf(s[T^). Then we have that 

| = &v,v' • (5.2) 


Clearly, 

E(||n^n,|| XH ( 7 )) =^( e ^\ u ) ||iW 0 >n„||, 


(5.3) 
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where E denotes the expectation value with respect to the measure 
deduce that 

£E(||n^">n„[|x a ( ?) ) <oc, 

n>0 


Using again (5.1), we 


(5.4) 


and hence that ||II !/ p( n )n ;/ ||x„(„) converges to zero, /Lt^-a.e.. For every element £ £ 3„/, there is 

“j,/ — 

no £ N such that £ £ for all n > no, and we conclude that \\H^ p^H ^||xs , converges to zero, 
/z^-a.e.. 


5.1.2 Completion of the proof of Theorem 3.3.1 


As (j 5 is finite and the relative entropy function I Pv is lower-semicontinuous, we can choose the 
sets { c ^v}ve’B. ao (also used in the previous paragraph) such that they satisfy: 


(a) n <% v , = 0 , 

(b) For every v ^ v' : I Pv ( p ) > min v ^ v , I Pv (p„>) -5/2 Vp £ 


where 0 < <5 < min I Pv {p u ' )- hr this particular case, it is convenient to identify „#( as ) with 
the compact convex subset of Rl <Ts l made of all vectors p = (pi,...,P| CTS |) such that pt > 0 and 


Y/iPi = 1 - We define again the sets Si," as the inverse image of the sets under the empirical 
Since 3 r 


measure yL £<r f^S^. Since 5^ is finite, we deduce from Theorem 
a £ (0,1) and C > 0 such that, for every v ^ z/, 


3.2.1 


that there are constants 


fc(5^H</a(4 n) )>)<cv\ 


Using (5.51 we obtain 


E 


(x s M J2 Tr(n< Ca n . 


v' 


Eq. (|5.6|) and some easy calculations imply further that 

e **<») TT(n„ P W)n„p^n„ - n„p( 0 ) n„ 


< CaP 


(5.5) 


(5.6) 


(5.7) 


Following similar procedures as above with x^M 
\\Hvp( n ' ) Il u \\x w (.'ni, using the fact that liminfS^"^ 


Tr(n u pWn v )n v p( n )n v 

= 3„, we find that 


n„p(°>n„ 


instead of 


lim xe.v 

n—> oo 


TX{Jl v p^J\ v )n v p^Yi v 


iW 0) n v 


= o, 


(5.8) 


/^-almost surely. 

5.2 Proofs of results in Section 0] 

In this part of the appendix we prove Lemma |4.2.1| For ease of comprehension, the proof is 
subdivided into a series of lemmas. The inequalities of Lemma |4.2.1| are proven in paragraphs 
[o~2~2l and PT~2~TT1 

5.2.1 Auxiliary lemmas 

Our first result concerns the speed of “purification” in non-demolition measurements. We recall 
that the functionals 

n 

Mu>(£l,-,fnM = IU(&l l/ ) and = ^2 Ti(IL v p i0) )ji u £ n \v) (5.9) 

i—1 t'ESoo 
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generate probability measures on the space (5, E) that we denote by the same symbols, and that 
the probability distributions p v = and TV = X^ g(rs p(£ W)^ on a S are assumed 

to be distinct unless v = v'. Furthermore it is assumed that p(£\u) 0 for all £, v, see Item (ii) of 
Assumption |4.1.1[ and that the sets 5^ and as are finite. 

Lemma 5.2.1. Let 6 £ (0, min„^„/ I Pv (tv)). There are constants C > 0 and a £ (0,1) (inde¬ 
pendent of p^), and a family of sets C S (independent of p^) such that the following holds 
true: 

(i) Hi n) nS^ = 0 and |/i w (s[ / ? ) \v)\ < Ca n if v ^ v'. 

(ii) For all v £ 


|/v(^” } ) - Tr(n,p(°))| < Ca n , 1 - fl u { (J 4 "' } ) < C\ 

(in) If ( £ si"-* and ifTr(p( 0 ^Il„) > 0, f/ien 


VoK- 


(5.10) 


E Tr(n^pW(C (n) )) < E 


Tr(n^/p(°)) - n(minI Ply ( Pl/ ,)-6) 


^ Tr(n,p(0)) 

v'v'r^v 


(5.11) 


Proof. The assumption that p((\v) 0 for all £, v implies that the sets a i/ v constructed in Para¬ 

graph |5.1.2| can be chosen such that the further property 

(c) For every 1 /, I p „(p) < f, Vp£^,, 


is satisfied (as well as (a) and (b) of Paragraph |5.1.2 1, because the relati ve entr opy functions are 
continuous in that case. The family of sets (s(," ) is defined as in Section 5.1.2 


/U^V) </U(4E>) < cv\ 

if v ^ v' and hence Item (i) is clear. Furthermore, 

£u,(5< n) ) - Tr(lW 0) )| = | E ^(W 0) )/U4"V) - Tr(II„/? (0) ) 


and still satisfies 
(5.12) 


< Ca n 


(5.13) 


where we have used that Y)u' T^IIj,/// 0 )) = 1. This implies the first equation in Item (ii). The 
second equation of Item (ii) follows easily: 


!-£<,( u ^ n) )H 5 Tr(n ^ (0)) ^ 5 ^( s - n) )l- c|E 

PGSoo PGSqo PGHqo 


(5.14) 


where we use (5.131. We finally prove Item (iii). We recall that the empirical frequencies are given 

by 

(5.15) 


/? } (i) =-#{.7 £{!,••• ,«} | &=£'}. 


We have that 


exp e = n p(?w) nf ^ ( - n 4 ' n) (o n/ ^ ,) 


(i) 


{'Go’s 


= /Uf‘“V) n 4A 


(5.16) 


25 



















and v' ^ v, 


Using (5.161, we deduce that for £ £ 
Tr ( n i / / / 5 ( n ) ( f (n) )) =^ 11 —IZ 


< 


Tr m „(».! < iKn-. P <°)) 

M£ (n> ) /UZZ) Tr(n y p(0)) 

Tr(n,^(o)) ) ex p( n/ p-( Z Z’(O^eO ~ n V( Z /e ( " } (£)*€'))• 


e'eirs 


As £ £ s£, n \ it follows from the property (c) that 

Trfn ,pW(fW)) < T r (n^p(°)) 
\ v p ’) s mflOD 


(5.17) 


□ 


Then we compare the density matrices p ^ and and their associated measures to density 
matrices obeying a non-demolition evolution. Given a density matrix p on JtFp and a finite length 
protocol £ <jg^, we introduce the random variable 

p(i (r) \R (k) ) ■■= ; r) o ... o $g 1} [p], (5.18) 

on <Tg r , where p : £ := (rfa, ..., rjk, £i,...,£;). We stress that p{P^\p^) is not normalized, and 
hence it is not a density matrix. 

Lemma 5.2.2. Let r, k £ N. For any set T C C7g r and for any £ erg^, we have that 

l|p(i (r) |?? (fe) )-^o...o$ 6 [p]|| < -^p—d-r- 1 Tr(p(p%W)), (5.19) 


Z |r? (fc) ))-TV(%. 

U r) er 


o ... o $ 


«1 


[p])) 


< 


dl ,- r -l 


d - 1 -1 


(5.20) 


where the first inequality holds for all £ <jg r . 

Proof: Assumption 4.1. l[ i) and a standard telescopic estimate imply that for any trace class 
operators pi,p 2 , 

ll*$>Pi - hM\ < llZlI (dilbill + IIPi ^ Pall) ■ (5.21) 

Iterating this inequality, we get that for any density matrix p, 


l|p(i (r) l2 (fc) )-%0"-o$ Cl [p]|| < Trm (r) \R (k) )) 


(5.22) 


where we have used that ||4>£p|| = Tr($^p) for all density matrices p , and that 

II^WlI • • ■ ||% + .- m+1 l|Tr($£ +r - m) o ... o $J+ 1} [p]) < dT m Tr(p(ZV fc) )), (5.23) 


which follows directly from repeated use of (4.6 1 and (5.211. To prove (5.201, it is sufficient to use 
that 

]TV(p(i (r V fe} )) - Tr(l>^ o ... o %[ p})\ < ||p(| M |p (fc) ) - $ £ „ o • ■ ■ o $ ?1 [p]|| 

and to sum over £ T. □ 
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5.2.2 Proof of the inequalities (4.12l and (4.14l of Lemma 4.2.1 


We use the family of sets E„^ C S defined in Lemmata 5.2.1 We introduce the family of sets 

~(k,k+r) — xfc cf(r) r w 

The projection on the first l entries of an element £ is denoted by piy: pr ; (£) = ^ l \ It clearly 
holds that 

-ETr(H v p^) 

= E( E r ^ 1 ’(4 | j I . o • • • o [p^ k \r/- k ^)]) — Ti(Il v p^ k \^ k ^))^p, u {r/' k ^). 

y( fc ) ^ r '>Gpr T .(zl r) ) 


Using Eq. (5.10) of Lemma 5.2.1 we deduce that 

|/i ta ,(3< fc - fc+r) ) - ETr(n„p<*>)| < Ca r . 
Now we notice that for sufficiently small 5, 


(5.24) 


^(fc.fc+r) r n _ i 

“i/ L ^y(fc,fc+r) = „ ana 


- 1 r yY ( < k ^ + r ) =1/ ^ V 1 "" 1 !/ 


c (e* m Y („' * 0, 


as follows from the construction of the sets see Section 5.2.1 From Eq. (5.121 using 

(5.9), arguing as above, we obtain that 

flu (5>(».»+,., = „ \ 3p +r) ) < CaU (5.25) 

Eqs. (5.24)-(5.25) imply (4.121. We now prove (4.141. By definition of i^ k ’ k+r ^(v) and from (5.91, 
we deduce that 

~e^ k +r)(y)= Y, Tr((l-n v )p {k+r \i (k+r) )) fi u (i (k+r) ) 

|6 = ^(fc,fc+r-) = „ 

= E Y Tr (W° } ) M^V)- 

v' C^‘ = ‘^ / y : (fc,fc + r) =u 

Eq. (4.14) follows by using that S y(k,k+r) =y C (S^’ fc+r - ) ) c if v ^ v' and Eq. (5.121. 

5.2.3 Proof of the inequalities (4.131 and ( |4.15| ) of Lemma 4.2.1| 

with ?y( fe ) = and p = pU) (^U)^ Furthermore, 


5.2.2 


We use the same notations as in Lemma 
we use the notation £( fe > fc + r ) = (£ fc+1) ...,£ fe+r ). We have that 

e (M+r)( l/ ) = Y Tr((l-n,)p( fe+r )(^ (fe+r) )) pU{ {k+r) ) 

“ J/ ( k, k +r) _ v 

Y Tr((i-n I/ )p(| (fe ’ fe+r) || (fe) )) pUi (k) ) 

l^ 3 Jf(.k,k+r ) =v 

Lemma [5 .2.2 1 implies that 

\\m ik ’ k+r) \i {k) ) - h k+ r ° -°E + > (fc E (fc) )]|| < Km ik ’ k+r) \i ik) ))’ ( 5 - 26 ) 

and hence we deduce that 

e {k ' k+r \v) < Y Tr ((! - n,)^ fe+r o • • • o $ £fc+1 [p< fc >(f<‘>)]) pUi (k) ) + ^ T rf- r ” 1 


i 6 =^(fc,fe + r) = 


fc+r 


< E E ^(i (fc) ) Tv(n^p (fe) (e (fe) )) II K&K) + iwErrf 


—i—1 


V'^V^E jp(k,k+r) =v 


i=k -\-1 


and (|4.15 l follows by using (5.121 and the fact that S ^ (fcife+r)=i/ C (S^’ fe+r )) c if i/ ^ v'. Eq. 
(4.131 follows from similar calculations and is left to the reader. 
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5.3 Some properties of completely positive maps 

In this section we prove the statement left out in Section |3.1| We consider a family of completely 
positive maps acting on a finite dimensional space such that = X^ eo . s < t ) *f is 

a unital map. We call a collection of operators T a satisfying 

^(as)[X] =J2KXT a 

a£l 

a Kraus decomposition of $*(<rs). Note that such a decomposition can always be taken as a sum 
of decompositions of over £, i.e. the index set I can be written as a union of sets I% with 
£ £ as and 

= £*.£[*] = ££ T* a XT a . 

The following preparatory lemma is well-known (see e.g. ED- 

Lemma 5.3.1. Suppose that <5>(crs) has a faithful stationary state p. Then kcr(<E>*(<7 s ) — 1) = 
{r«}U/ for any Kraus representation T a . 

Proof: We recall that $(o\s) = (<h*(ers))* and that p is a stationary state if <E>((7s)[p] = p and 
a faithful state if Tr{pA* A) = 0 implies A* A = 0. We write = $*(<r s ). 

If A £ ker(<h* — 1) then so does A*. Then Kadison inequality $*[AA*] > [A]$*[A*] = AA* 

holds and the difference between the LHS and the RHS can be explicitly expressed as 

<S>4AA*}~ AA* = |[-A*,r Q ]| 2 . 


Applying the stationary state, Tr(p-), on both sides we get Tr(p^) Q |[A*,T ct ]| 2 ) = 0. Since p is a 
faithful state the statement follows. □ 

The following lemma gives a description of a non-demolition family $*£. 

Lemma 5.3.2. Suppose that 

o $*£/ = o (5.27) 

for all (f £ as and that the map <h(<rs) has a faithful stationary state. Then $* £ '[11]] = 0 

and the corresponding joint spectral decomposition of the identity 11 = n„ has the property that 

<fr*f[n„] = p(f |z/)n„ for some p(£\v) > 0. 


Moreover if v ^ v' then there exists f such that p(£\v) ^ p{f\v'). 


Proof: Since $*(crs) = X] £e<7s ( |5.27| ) implies that $*(as) o $*£ 

<&*(crs) o $* £ [11] = <l?* £ [ll]. 


o <h*(crs) and hence 


It then follows from the previous lemma that [$*^[1], Tq,] = 0 for any Kraus decomposition T a of 
<I>*((Ts). We then have 


$**' 0 $ *€pL] = E r a^W r a= E 

Repeating the same with £ and £' exchanged we conclude that and commute. Let 

now 11 = II„ be the joint spectral decomposition of the commuting family {<!>*£[U]}£ eo - s and 
suppose that 

<j>* ? [i] = Ep(£K)n,'. 
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By hitting the equation by II„ from the right and using the commutativity of the operators r a 
with n„ we get 

= p(£\v)n v . 

To finish the proof notice that p(£ \u) is non-negative because $*£ is a positive map. Furthermore, 
if p(£ \v) = p(£|i/) for all £ then the projections !!„ and II„' would not appear separately in the 
spectral decomposition, but rather their sum II„ + n„/ would appear. □ 
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